Pressure-induced magnetic transition and volume collapse in FeAs superconductors: 

An orbital-selective Mott scenario 
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Motivated by pressure experiments on CaFe2As2, we propose a scenario based on local-moment 
physics to explain the simultaneous disappearance of magnetism, reduction of the unit cell vol- 
ume, and decrease in resistivity. In this scenario, the low-pressure magnetic phase derives from Fe 
moments, which become screened in the paramagnetic high-pressure phase. The quantum phase 
transition can be described as an orbital-selective Mott transition, which is rendered first order 
by coupling to the lattice, in analogy to a Kondo volume collapse. Spin-fluctuation driven super- 
conductivity competes with antiferromagnetism and may be stabilized at low temperatures in the 
high-pressure phase. The ideas are illustrated by a suitable mean-field analysis of an Anderson 
lattice model. 
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I. INTRODUCTION 

The discovery of superconductivity in the iron 
pnictides, 1-3 with the maximum reported transition tem- 
perature over 50 K 4 has spurred intense activities in 
both experiment and theory. While the initial re- 
search has mainly concentrated on compounds with 
PA/nmm structure, e.g. LaFeAsOi-^F^,, superconduc- 
tivity was later also found in oxygen-free materials with 
I4/mmm structure, 5 with examples being CaFe2As2 and 
BaFe2As2- According to their chemical composition, the 
two groups have been dubbed 1111 and 122 compounds, 
respectively. In both cases, the structure consists of FeAs 
layers, which are believed to be responsible for the low- 
energy electronic properties. 

The temperature-doping phase diagram of 
LaFeAsOi-zF-r has been mapped out in some detail. 6 At 
x = the material undergoes a structural transition into 
an orthorhombic phase upon cooling, closely followed by 
a magnetic transition. Neutron diffraction shows that 
the order is of layered antiferromagnetic type, with an 
in-plane ordering wavevector (tt,0)J With increasing x, 
magnetism disappears abruptly at x ~ 4.5% and gives 
way to superconductivity, with T c having little doping 
dependence up to x — 0.2. 

This phase diagram suggests interesting similarities to 
the high-T c cuprates, where a doping-driven transition 
from an antiferromagnet to a superconductor is found as 
well. However, a number of differences between the FeAs 
compounds and the cuprates are apparent, (i) The "un- 
doped" FeAs materials are not Mott (or charge-transfer) 
insulators, but may be characterized as bad metals, (ii) 
The FeAs magnetism is not fully consistent with a local- 
moment picture, as the size of the ordered moment is 
unexpectedly small, in the range 0.25 ^_b 8 to 0.36 hb, 7 
whereas a local moment on Fe is expected to have at 
least 2 fis ■ Although magnetic frustration has been iden- 
tified as a relevant ingredient, 9 this alone is not enough 
to explain the experimental data, (iii) In the doping 
regime where superconductivity occurs at low tempera- 



tures in the FeAs materials, the signatures of magnetism 
appear to be weak, although the magnetic susceptibil- 
ity is larger than in a conventional metal, 10 and a mag- 
netic "resonance" mode has been reported below T c in 
Ba(Fe,Co) 2 As2. n 

From first-principles calculations, 12 the most impor- 
tant bands in the iron arsenides are of Fe 3d charac- 
ter, with sizeable admixtures of As p states. Impor- 
tantly, all Fe 3d orbitals appear to contribute to the low- 
energy properties. In contrast, in the cuprates only the 
Cu id x 2_ y 2 orbital (together with O 2^ x , v ) contributes. 
Thus, orbital degeneracy and strong hybridization likely 
play an important role in the iron arsenides. (It should be 
noted that the reported results of first-principles calcula- 
tions are partially contradictory, in particular regarding 
the magnetic properties, see Ref. 13.) 

Pressure experiments on CaFe2As2 have provided in- 
teresting additional information. 14-19 Starting from the 
undoped magnetically ordered compound in the or- 
thorhombic phase, the application of hydrostatic pres- 
sure induces a first-order transition at around 0.4 GPa 
into a (compressed) tetragonal phase, where magnetism 
disappears. Remarkably, there is also a pressure-driven 
structural transition at elevated temperatures, from the 
ambient-pressure tetragonal phase into the compressed 
high-pressure tetragonal phase. This transition is ac- 
companied by a drop in the resistivity, i.e., the high- 
pressure phase appears to be a better conductor. 15 Band- 
structure calculations suggest that the disappearance 
of magnetism is intimately connected to the pressure- 
induced volume collapse. 14 ' 16 ' 20 The high-pressure phase 
was reported to be superconducting 14,15 with a T c of 
about 12 K over a significant range of pressures, how- 
ever, subsequent experiments 17,19 with optimized hydro- 
static pressure conditions did not detect superconductiv- 
ity down to 4.2 K, suggesting that the earlier reports were 
related to strain-induced superconductivity. Neverthe- 
less, all measurements agree on the volume collapse and 
the concomitant loss of magnetism and drop of resistiv- 
ity. In addition, superconductivity under high pressures 
has also been reported 21 in SrFe2As2 and BaFe2As2- 
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The purpose of this paper is to propose a phenomcno- 
logical local-moment-based scenario for the physics of the 
FeAs compounds, with the primary goal of explaining the 
pressure experiments on CaFc 2 As 2 . 14-16 ' 19 Our scenario 
is based on the assumption that bands with more itin- 
erant electrons co-exist and interact with more localized 
(i.e. moderately to strongly correlated) ones, which ren- 
ders the problem similar to Kondo or Anderson models 
for heavy- fermion metals. (We note that related ideas on 
localized and itinerant electrons in iron arsenides have 
been put forward in a few recent theory papers 22-25 which 
appeared while this work was being completed.) Our 
analysis will draw analogies to recent developments in 
heavy-fermion physics, particular, we shall associate the 
pressure-induced transition in CaFe2As2 with a variant of 
the Kondo volume-collapse transition, driven by a large 
increase of the hybridization between the itinerant and 
localized electron bands, which in turn quenches the mag- 
netism. 

The remainder of the paper is organized as follows: In 
Sec. II we sketch our general ideas and point out sim- 
ilarities and differences between the iron arsenides and 
classical heavy-fermion materials. We then illustrate the 
proposal by a simple mean-field calculation in Sec. Ill, 
which we believe captures important parts of the rele- 
vant physics of the CaFe2As2 pressure experiments. A 
discussion of implications will wrap up the paper. 



II. LOCAL MOMENTS IN A CORRELATED 
ANDERSON LATTICE 

Our proposal is based on the assertion that electronic 
correlations in the iron arsenides are sizeable. 9 ' 22,26 ' 27 
Then, the magnetism is not purely of weak-coupling type, 
but instead local-moment physics is relevant. 28 This view 
appears consistent with the results of neutron scattering 
experiments, which find spin-wave excitations in the an- 
tiferromagnetic ground state of CaFe2As 2 , that are well- 
described by an anisotropic 3d Heisenberg model, at least 
for long wavelengths. 29 

To explain the fact that the antiferromagnetic phase 
of FeAs compounds is neither an insulator nor a good 
metal, we invoke the existence of (at least) two types 
of electrons, one more localized and one more itiner- 
ant species. In the spirit of a Kondo or Anderson lat- 
tice model, we shall adopt here the language of "local- 
moment" and "conduction" electrons. In the antiferro- 
magnetic phase, the local-moment electrons carry the 
magnetism, and the residual interaction between local- 
moment and itinerant electrons provides sizeable scat- 
tering, leading to both bad-metal behavior and reduced 
moment amplitudes. 28 ' 30 In the paramagnetic phase, the 
local-moment electrons are strongly hybridized with the 
conduction electrons (i.e. the moments are Kondo- 
screened), leading to good metallic behavior. Provided 
that spin fluctuations are still sizeable, spin-fluctuation 
mediated superconductivity can arise at low tempera- 



tures. 

How does this phenomenological picture tie in with 
microscopic considerations? The most plausible scenario 
is that both types of electrons are primarily of Fe 3d 
character, with strong admixtures of As 2p. As detailed 
in Ref. 22, the interplay of p-d hybridization, spin-orbit 
coupling, and crystal-field splitting can lead to the two 
highest occupied FeAs levels being filled with one elec- 
tron, but with different correlation strengths. The elec- 
trons in the lower level are more localized, providing a 
natural basis for our phenomenological approach. 31 

Assuming the applicability of an Anderson lattice pic- 
ture, let us further develop our ideas, concentrating on 
the transition between the antiferromagnetic and param- 
agnetic phases. In fact, many qualitative results can be 
borrowed from recent theoretical work on heavy-fermion 
materials. There, the zero-temperature transition be- 
tween a paramagnetic heavy Fermi liquid and an an- 
tiferromagnetic metal can be a standard spin-density- 
wave transition of a Fermi liquid, described by a Landau- 
Ginzburg- Wilson theory of Hcrtz-Millis type, or involve 
the breakdown of Kondo screening. 32 For the latter sce- 
nario, a number of theoretical descriptions have been put 
forward. 33-37 Within the Kondo-breakdown scenario of 
Scnthil et al., 36 the coupling to lattice degrees of free- 
dom has been investigated recently. 38 If the electron- 
lattice coupling is sufficiently strong, then the Kondo- 
breakdown transition is rendered first order and accom- 
panied by a isostructural volume change, which is a zero- 
temperature variant 39 of the classical Kondo volume col- 
lapse transition. 41-44 This first-order transition extends 
to finite temperatures and masks an otherwise existing 
magnetic quantum critical point. (A first-order transi- 
tion removes the sharp distinction between the Kondo- 
breakdown and spin-density-wave transition scenarios.) 

In a situation with intermediate correlations, where 
valence fluctuations are not fully quenched, it is useful 
to think about the Kondo breakdown transition as an 
orbital-selective Mott transition: 37 In the non-magnetic 
Fermi liquid phase, the local-moment electrons become 
itinerant and strongly hybridized with the conduction 
electrons, whereas they undergo Mott localization in the 
magnetic phase. In fact, such orbital-selective Mott tran- 
sition have been studied extensively in two-band Hub- 
bard models. 45 Based on the cited works, we propose 
that the pressure-driven transition in CaFe2As 2 can be 
described as orbital-selective Mott transition, which be- 
comes strongly first order due to electron-lattice cou- 
pling. The low-pressure phase displays partial Mott lo- 
calization, leading to local-moment magnetism, whereas 
in the high-pressure phase strong hybridization quenches 
the moments and leads to a conventional paramagnetic 
Fermi liquid. The pressure-driven iso-structural lattice 
transition at elevated temperatures is then the analogue 
of a Kondo volume collapse. 

It is important to point out that our scenario is not 
in contradiction with itinerant spin-density- wave descrip- 
tions of magnetism in the FeAs compounds. In the 
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Kondo-lattice context, itinerant and local-moment mag- 
netism can be adiabatically connected (apart from transi- 
tions involving changes of the Fermi surface topology). 30 
Therefore, a system at intermediate coupling can in prin- 
ciple be described using both itinerant and localized elec- 
tron concepts. Here, we find it advantageous to employ 
a strong-coupling language in order to highlight the na- 
ture of the volume collapse transition. Independent of 
the language, the driving force of the volume collapse 
is the energy gain due to an increasing in the effective 
hybridization between the bands. 

Let us note that there are a few important differences 
between our envisioned scenario for the FeAs materials 
and the phenomenology of heavy- fermions compounds, 
(i) The local-moment electrons in the iron arsenides are 
probably far from the Kondo limit. Then valence fluctu- 
ations are sizeable, and the picture of Kondo screening 
does not literally apply. The coherence temperature on 
the high-pressure side of the transition is not small, and 
the quasiparticles are not very heavy, in contrast to that 
of typical heavy fermions. (ii) The character of the (it, 0) 
magnetic order is important to understand the details of 
the phase diagram. 

The magnetism deserves a few further comments: In 
order to explain the magnetic order at wavevector (n, 0), 
both local-moment and itinerant scenarios have been in- 
voked (which are not mutually exclusive. 28 ' 30 ) In a local- 
moment picture, the exchange interactions J\ and Ji on 
the square lattice of Fe atoms (between nearest and next- 
nearest neighbors, respectively) have been deduced to be 
both antiferromagnetic with Ji < J2. 9 ' 12 In this regime, 
the 3\-Ji square-lattice Heisenberg antiferromagnet is 
known to have a "layered" antiferromagnetic ground 
state with (n, 0) order. 46 ' 47 In contrast, in an itinerant 
picture the (tv,0) order arises from nearly nested Fermi 
surface pieces. Independent of whether the magnetism 
is better described in an itinerant or localized picture, 
magnetic ordering at (tt, 0) in a originally tetragonal en- 
vironment breaks the 90° lattice rotation symmetry. This 
induces an orthorhombic distortion inside the antiferro- 
magnetic phase. However, the lattice rotation symmetry 
may also be broken at a higher temperature than the spin 
symmetry, in which case the orthorhombic distortion oc- 
curs before the magnetic order. 48,49 This is indeed what 
happens experimentally, i.e., magnetic fluctuations are 
likely the driving force of the structural phase transition. 

Last not least, we discuss the possible emergence of 
superconductivity. For the Kondo-brcakdown scenario, 
magnetically mediated pairing has been argued to be a 
generic instability of the Fermi liquid, 35 with the max- 
imum T c near the continuous Kondo-breakdown transi- 
tion. In the present situation of intermediate correla- 
tions, the absolute value of T c will depend strongly on 
microscopic details, e.g., of the band structure. More- 
over, if the lattice coupling renders the transition strongly 
first order, the system will effectively "jump" over the pa- 
rameter region with large T c . Thus, while both pressure- 
induced volume collapse and loss of magnetism are inte- 



gral and robust parts of the proposed scenario, supercon- 
ductivity with a sizeable T c is more fragile. 

III. MEAN-FIELD THEORY 

We now illustrate the ideas described in the last section 
by a simple model calculation. We shall refrain from us- 
ing a realistic band structure with five or more bands, but 
instead employ a two-band Anderson lattice model which 
is sufficient to capture most of the qualitative physics. 

A. Anderson-Heisenberg lattice model 

The starting point of our analysis is an Anderson lat- 
tice model, describing delocalized conduction (c) elec- 
trons on a lattice which hybridize with correlated and 
more localized / electrons on the same lattice. (Despite 
the labels c and /, both bands may have primarily Fe 
3d character.) To simplify the approximate treatment 
of magnetism, we supplement our model by an explicit 
Heisenberg exchange interaction between the local mo- 
ments. The resulting Hamiltonian reads 

ka icr i 

kia 

+ J2 J *(h3)S t -S 3 (1) 

ij 

in standard notation. The first term describes conduc- 
tion electrons with band filling n c , £° (U) is the bare / 
electron energy (Coulomb repulsion), Vk the hybridiza- 
tion matrix element, and M the number of unit cells. 
Si = fiu^aa'fia' /2 is the operator for the spin moment of 
the / electrons on site i. Although the system is three- 
dimensional, we shall neglect the electronic coupling be- 
tween the layers for simplicity, hence the model (1) will 
be treated on a 2d square lattice. Further we shall assume 
Vk = V°, and take the Heisenberg interaction Jn(i,j) to 
be non-zero for nearest and next-nearest neighbors, with 
values J 1 < J . Note that the Kondo limit will not be 
taken. 



B. Lattice distortions and electron— lattice coupling 

Elastic energy changes of the lattice under application 
of hydrostatic pressure depend on the lattice distortion. 
For tetragonal lattice symmetry, the most general form 
of the elastic contribution to the free enthalpy is 40 

Gi at (e x ,e y ,e z ) = l;v c 3 el + ^v ci(e 2 x + e 2 y ) + c 12 e x e y 
+ v Q c 13 e z (e x + e y ) + pv (e x + e y + e z ){2) 
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Here, the dimensionless e x ,y.z are the diagonal entries of 
the strain tensor, i.e., the relative changes of the lattice 
parameters of the tetragonal unit cell. The elastic con- 
stants Cj, Cij, and the reference volume depend on the 
material at hand. Below, we shall choose values for the 
d, dj such that the experimentally observed lattice dis- 
tortions of CaFe2As2 (Ref. 14) are approximately repro- 
duced. Importantly, the largest pressure-induced change 
is the collapse of e z , whereas e x , y even increase slightly 
at the collapse. 

We now construct a model for the electron-lattice cou- 
pling, combining the results from different approaches 
with simple theoretical arguments. We start with the 
strain dependence of the magnetic couplings J\ and 
J2 which are responsible for the ambient-pressure mag- 
netism, and will also be relevant for spin-fluctuation- 
mediated superconductivity. In experiment, antiferro- 
magnetic ordering at wavevector (it, 0) is energetically 
stabilized by an orthorhombic distortion with e x > e y . 
Additional evidence is the strong dependence of J\ 
and J 2 on the strain e z , as shown in band structure 
calculations. 50 In addition, it has been shown 56 that 
changes in the Fe-As bond angle stabilize antifcrromag- 
nctic ordering at wavevector (tt, 0), providing an addi- 
tional hint towards the importance of coupling all strain 
directions to the exchange couplings. We therefore pa- 
rameterize 



is 



Jix = ■/i(l+7i"e*+7?e*), 
Jiy = ^i°(l + 7i%+7i%), 
J 2 = 4 [l + li{e x + e y )+ 1 z 2 e z ] 



(3) 



From neutron scattering data and LDA calculations for 
the orthorhombic phase of CaFe2As2 29 a rough esti- 
mate of 7^ can be obtained. For S = 1/2, we plug 
the LDA values J\ x — 82 meV and J\ y = 20 mcV 
from Ref. 29 into the parameterization of Eq. (3) and 
set e x — e y = 0.01, as observed in the orthorhombic 
phase of CaFe2As2- 14 This calculation yields the estimate 
j 1 w 70. The size and sign of 7^ might well account for 
the observed orthorhombic distortion and the (n, 0) or- 
dering vector, as argued in Ref. 12. From band-structure 
calculations for LiFeAs 50 it is suggested that and 7| 
are numbers of O(10), whereas their signs turn out to 
depend on pressure. To account for the vanishing super- 
conductivity at high pressures, we chose their sign to be 
positive, since 7* and 7I dominate the changes in J\ and 
J 2 at high pressures. 51 Finally, 7^ is of subleading influ- 
ence within our calculation as long as it is of O(10) or 
smaller, which is suggested by the order of magnitude of 
7f and 7|. 

A crucial ingredient for our model is the strain depen- 
dence of the hybridization, which we assume to be the 
dominant mechanism to drive the volume collapse, as ob- 
served in the Kondo volume collapse in Cerium and other 
materials. In general, the mechanism for the Kondo vol- 
ume collapse transition is a gain in hybridization energy 
via a structural distortion. A plausible parameterization 



V = V° [l+^(e x + e y ) +1 z e z 



(4) 



As our model is to be understood as an effective model, V 
is related to actual band structure parameters in a likely 
complicated fashion, and information about its strain de- 
pendence is not available. From the Kondo volume col- 
lapse model it is known 39,42 that 7- 1 , j z should be cho- 
sen of O(l) to reproduce volume collapses of O(10%). 
Given the experimental behavior of the CaFe2As2 lattice 
constants, 14 7 Z < is crucial. We shall also take r y- L > 0, 
but this is subdominant. 

Finally, we neglect any dependence of the c dispersion 
£k and the / energy eo on lattice strain. While such a de- 
pendence certainly exists, it will not qualitatively change 
our picture, which is dominated by the effects parame- 
terized in equations (3) and (4); moreover, a pressure de- 
pendence of the bandwidth can be absorbed in a pressure 
dependence of the reference energy scale of our calcula- 
tion. 



C. Fermionic mean- field theory 

The model (1) will be solved using a standard mean- 
field theory, with a fermionic representation of the local 
moments and slave bosons to deal with strong local repul- 
sions. Theories of this type have been extensively used 
in the study of the Kondo and Anderson lattice and also 
allows to deal with the magnetic exchange term in Eq. 

Q\ 36,37,52,53 

In the limit of infinite Coulomb repulsion U, the phys- 
ical electron can be represented by a spinless boson j-j 
and an auxiliary fermion f ia , f icr — r\j iai together with 
the constraint 



+ f f f = 1. 



(5) 



The auxiliary bosons rj will be treated on the mean- 
field level. The Heisenberg part of the Hamiltonian, rep- 
resented in pseudo-fermions ji a , has to be decoupled. 
Guided by the presence of potential magnetic and su- 
perconducting instabilities, we introduce mean fields ac- 
cording to 



Mr = la 



' fra') t 



A i] - -(fiffjl - filfjt)- 



(6) 



The decoupling of the quartic interaction takes the form 
(additional constants omitted): 



1 

4 



Si ' Sj — X 

+ (I-*) 



-A l3 (fyi-fj l f],) + h. c . 



t jt 



•(7) 



We have introduced an additional decoupling parameter 
x e (0, 1), where x = corresponds to an Sp(iV) large- 
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limit 35 and x = 1/2 to unrestricted Hartree-Fock. Physi- 
cally, x regulates the balance between ordered magnetism 
and superconductivity, and we shall choose x = 0.3 and 
0.4 below. 

It remains to specify the spatial dependence of the 
mean-field parameters. The obvious parameterization 
for the magnetization M r is M r = m s exp(iQ • r) with 
Q = (7r,0). 54 The complex pairing fields Aj 3 - live on the 
nearest-neighbor and next-nearest-neighbor bonds of the 
square lattice. In the general orthorhombic case, we fo- 
cus on saddle points with A x on horizontal bonds, — A y 
on vertical bonds, and ±A xy on diagonal bonds, such 
that the pairing term in fc-space reads Ag = A x cos k x — 
A y cos k y + 2A xy sin k x sin k y . For a tetragonal lattice, 
A x = A y , and the pairing is a mixture of d x 2_ y 2 and d xy 
symmetric terms. 54 Finally, the slave boson is chosen to 
be uniform, rj = ro. 

The mean-field amplitudes are obtained from the 
saddle-point equations 



A., 
A„ 



COS ktr 



4M 



+QT 



fklfk+Ql 



1-ri 



(8) 



where A denotes the complex conjugate of A. 



D. Phases and electronic phase diagram 

The electronic mean-field theory specified above can 
display the following phases (not all of which will appear 
for our choice of parameters) : 

• Decoupled, with ro = 0, m s = 0, A = 0, describ- 
ing a paramagnetic high-temperature regime with 
localized / electrons. 

• Fractionalizcd Fermi liquid (FL*,) with ro = 0, 
m s = 0, A / 0. This is a paramagnetic phase, 
where conduction electrons alone form a "small" 
Fermi surface (FS) and are decoupled from a frac- 
tionalized spin liquid of paired spinons. This phase 
was introduced in Ref. 35, but will not play a role 
here. 

• Local- moment antiferromagnet (AFM), with r = 
0, m s 7^ 0. Here, A may be zero or finite, the latter 
case reflecting residual spinon pairing. 
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FIG. 1: Electronic mean-field phase diagram in the 
temperature-hybridization plane, for fixed lattice parameters 
£x,y,z = 0. Top: Mean-field decoupling parameter x — 0.3. 
Bottom: x = 0.4. Thick (thin) lines denote first-order (con- 
tinuous) phase transitions. For small hybridization, the / 
electrons are localized, and magnetism dominates at low T, 
whereas large hybridization leads to itinerant / electrons and 
superconductivity. For details and parameters see text. 



• Paramagnetic Fermi liquid (FL), with r ^ 0, m s — 
0, A = with itinerant / electrons and a "large" 
Fermi surface. 

• Antiferromagnetic Fermi liquid, with ro ^ 0, m s ^ 
0, A = with itinerant / electrons. 

• Paramagnetic superconductor (SC), with ro 7^ 0, 
m s = 0, A 7^ 0, obtained from pairing in the large- 
FS Fermi liquid. 

• Antiferromagnetic superconductor, with ro =/= 0, 
m s / 0, A / 0. 

A zero-temperature transition from a phase with ro 7^ 
to ro = is associated with Mott localization of the 
/ electrons, i.e., an orbital-selective Mott transition. A 
continuous T = transition of this type will survive be- 
yond mean-field in the paramagnetic case, 35,36 whereas 
it becomes a crossover in the antiferromagnetic case. 30 

For the numerical calculations, serving as an illus- 
tration of our ideas from Sec. II, we choose parame- 
ters as follows: The c electron dispersion consists of 
nearest-neighbor hopping, Ek — — 2£(cos fc^+cos k y ), with 
t = 0.5 eV and filling n c = 0.8. The resulting band width 
is comparable to the energy range where the Fe DOS is 
sizeable. 55 The / level position is set slightly below the 
lower band edge to obtain moderate valence fluctuations, 
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ff = — 2.3eV. The properties of the local-moment anti- 
ferromagnetic phase are determined by the exchange cou- 
plings xJ\ and xJi- For a stable AFM ground state with 
wavevector (ir,0), we choose the ratio J^/^i = 1-5- We 
employ two parameter sets: x = 0.3, Jf = 200 mcV, Ji? — 
300 meV and x = 0.4, J? = 150 meV, J 2 ° = 225 meV. 
The resulting x J® and xJ% coincide with the theoretical 
values for BaFe2As2 from Ref. 57. For CaFe2As2, exper- 
imental and theoretical values have been determined in 
the orthorhombic phase in Ref. 29. For our model with 
5=1/2, those results yield Ji « (J lx + J ly )/2 ss 40meV 
and J2 50meV. 

The electronic phase diagram, obtained for fixed 
£ x,y,z — as function of temperature T and hybridization 
Vo, is shown in Fig. 1, for decoupling parameters x = 0.3 
and 0.4. We find magnetism and superconductivity to be 
mutually exclusive and separated by a first-order transi- 
tion: At this transition, both m s and tq jump, i.e., the 
system switches from a local-moment dominated antifer- 
romagnet to a fully itinerant superconductor. At low 
T, spinon pairing co-exists with local-moment antifer- 
romagnetism on the small- y side of the phase diagram. 
The thermal magnetic transition is very weakly first order 
within our accuracy, but we cannot exclude it to be con- 
tinuous. The superconductivity is of d xy character, i.e. 
driven by the exchange interaction J2. (A small id x 2_ y 2 
admixture develops at low temperatures for x — 0.3.) 
The finite-temperature transition between the decoupled 
regime and the FL will be smeared into a crossover by 
fluctuations beyond mean field; the other transitions are 
accompanied by physical symmetry breaking and survive. 

A variation of the decoupling parameter x within the 
range x € (0.25,0.45) (keeping xJ\ 2 fixed) mainly influ- 
ences the stability of the superconducting phase, Fig. 1. 
Increasing x disfavors pairing and further stabilizes an- 
tiferromagnetism. For small x < 0.25 a paramagnetic 
FL* phase is realized at small hybridization Vq, whereas 
for large x > 0.45 superconductivity disappears. For 
x > 0.5, antiferromagnetism appears even in the large-FS 
Fermi-liquid regime. (A related mean-field theory with 
spinon hopping (instead of pairing) generically displays 
an itinerant antiferromagnetic phase. 36 ) In general and 
beyond simple mean-field approximations, the supercon- 
ducting T c will sensitively depend on band structure de- 
tails and nesting conditions; a variation of x in our cal- 
culations mimics such changes. 



E. Phase diagram with electron— lattice coupling 

The central result of our mean-field study is the phase 
diagram in Fig. 2, which accounts for lattice distortions 
and external pressure. It has been obtained from min- 
imizing the free enthalpy G e \ + Gut, where G\ & t is in 
Eq. (2) and G e \ is the electronic contribution according to 
the Hamiltonian 7i (1), with the lattice dependence of all 
parameters as in Sec. Ill B and the mean-field approxima- 
tion as in Sec. Ill C. As experimental data for the elastic 
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FIG. 2: Mean-field phase diagram including electron-lattice 
coupling in the temperature-pressure plane. Top: Mean-field 
decoupling parameter x = 0.3. Bottom: x = 0.4. Thick 
(thin) lines denote first-order (continuous) phase transitions. 
The isostructural volume-collapse transition displays a critical 
end point at a temperature of roughly 200 meV. See text for 
details. 



constants of the 122 materials were not available to us, we 
choose elastic constants of c\ = 441 kBar, C3 = 198kBar, 
and C12 = C13 = 66 kBar. The employed electron-lattice 
couplings are 7 1 - = 3.1, -f z — —5.2, 7^ = 35, 7* = 1.0, 
J2 = 8.0, and 7I = 5.0, see Sec. IIIB. The values for Cj, 
Cij, 7" 1 , and r y z were adjusted to reproduce the experi- 
mentally observed lattice distortions. The couplings 7^, 
^2 and 7f had little influence on the lattice distortions. 
They were adjusted to reproduce the experimentally ob- 
served phase boundary of the superconducting phase. 

With this parameter choice, the main effect of increas- 
ing pressure is an increasing hybridization, while the 
magnetic exchange couplings decrease somewhat. Con- 
sequently, the topology of the phase diagram is similar 
to that of Fig. 1, however, with a few crucial differ- 
ences: (i) The orbital-selective Mott transition becomes 
strongly first order, at both low and high temperatures, 
(ii) The antiferromagnetic phase is accompanied by an 
orthorhombic lattice distortion, and the thermal phase 
transition is of first order, (iii) Pressure drives a volume 
collapse transition, which is tetragonal <-> tetragonal at 
elevated T and orthorhombic <-> tetragonal at low T. The 
pressure dependence of various microscopic parameters is 
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FIG. 3: Pressure dependence of several microscopic parame- 
ters for low T = 0.5 meV and decoupling parameter x — 0.3. 
a) Lattice parameters e x , y ,z- Below temperatures of 30meV, 
the e depend only weakly on temperature, therefore this plot 
is representative for most parts of the phase diagram, b) 
Magnetic exchange constants J\,2- c) Hybridization V. d) 
Pairing fields A. Note that A reflects superconductivity only 
in the Fermi-liquid regime at large pressure, whereas it only 
describes spinon pairing in the orbital-selective Mott phase at 
small pressure. 



features of the phase diagram. 

The quantitative values of the lattice distortions in 
Fig. 3a are close to the experimental ones. The pres- 
sure dependence of the magnetic exchange, Fig. 3b, is 
relatively strong and not confirmed experimentally, al- 
though the qualitative behavior is suggested by neutron 
scattering data and band structure calculations. In our 
calculation, the pressure dependence of J\$, is required 
to create a pressure dependence of both the Neel tem- 
perature and the superconducting T c comparable to ex- 
periment. (In principle, it is conceivable that fluctua- 
tion effects not captured by our treatment have a similar 
effect.) The hybridization V displays a large jump at 
the critical pressure, Fig. 3c - this is, within our theory, 
the mechanism which balances the increase in elastic en- 
ergy connected with the volume collapse. The pairing 
fields undergo several changes as function of pressure, 
see Fig. 3d for x — 0.3, which depend sensitively on the 
pressure dependence of the exchange couplings and other 
microscopic parameters. The dominant pairing is of d xy 
symmetry as above, both for the spinon pairing of local- 
ized / electrons in the magnetic small-pressure phase as 
well as for the superconductivity of itinerant / electrons 
at high pressure (again with a small id x 2_ y 2 admixture at 
low temperatures, which disappears for x = 0.4). Slight 
changes of Jnj J\ stabilize saddle points with other pair- 
ing symmetries, but leave the non-superconducting part 
of the phase diagram essentially unchanged. Finally, the 
/ occupation jumps from 1 to w 0.80 across the transition 
for all temperatures below 60meV (not shown). 



shown in Fig. 3 for low temperatures. 

Our results reproduce salient features of CaFe2As2 as 
reported in Refs. 14,19. In particular, the phase diagram 
Fig. 2 is very similar to Fig. 3c of Ref. 14. For x = 0.3, 
our superconducting transition temperature close to the 
critical pressure is of order 50 K, whereas the Neel tem- 
perature is somewhat below room temperature. While 
these values could be brought closer to the experimental 
ones of Ref. 14, we refrain from doing so, because inter- 
plane coupling and fluctuation effects, both absent from 
our treatment, will significantly modify the mean- field re- 
sult. For x = 0.4 the pairing fields are much smaller and 
the range of superconductivity is strongly reduced; for 
x = 0.45 (not shown) T c is below 0.1 meV, broadly con- 
sistent with the more recent experiments of Ref. 17. All 
other parameters behave essentially identically as com- 
pared to x = 0.3. Hence, we conclude that the volume 
collapse and orbital-selective Mott transition are robust 
features of our model, but, in contrast, superconductivity 
with a sizeable T c is not. 

The present mean-field approach is not able to distin- 
guish the nematic transition from the magnetic transi- 
tion, which are separated in experiment, but coincide in 
Fig. 2. This deficiency could be in principle repaired by 
including an additional nematic order parameter in the 
mean- field treatment, but would not change the other 



IV. DISCUSSION 

In this paper, we have proposed a theoretical sce- 
nario to rationalize the pressure-induced phase transi- 
tions in CaFe2As2- Underlying the scenario is the physics 
of the Anderson lattice, used for heavy-fermion metals: 
Strongly localized electrons hybridize with more itiner- 
ant electrons. Microscopically, both may be primarily of 
Fe 3d character, as discussed in Ref. 22 (although this is 
not required within our phenomenological approach) . At 
ambient pressure, the localized electrons order antiferro- 
magnetically in a collinear arrangement at low temper- 
atures, accompanied by an orthorhombic lattice distor- 
tion. The system is weakly metallic, due to the presence 
of the itinerant carriers with a "small" Fermi volume. 
Increasing pressure drives a transition towards a param- 
agnetic Fermi liquid, where the previously localized elec- 
trons become itinerant and non-magnetic. The coupling 
to the lattice degrees of freedom in CaFe2As2 renders this 
transition strongly first order - this is akin to a Kondo 
volume collapse (although the system is not in the Kondo 
regime, and valence fluctuations are sizeable). The first- 
order nature of the transition, already in the purely elec- 
tronic theory, implies a tendency towards phase separa- 
tion, which appears to be present experimentally. 58 

This set of ideas holds plausible explanations for (i) 
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the coincidence of volume-collapse and magnetic-non- 
magnetic transition (which is also borne out by first- 
principles approaches 14, 16,20 , (ii) the bad metallic behav- 
ior, and (iii) the reduced magnetic moment of the mag- 
netic phase. Note that (ii) and (iii) are not described by 
the mean-field theory, but arise from residual scattering 
between itinerant and localized electrons in the orbital- 
selective Mott regime as small pressure. Moreover, the 
abrupt disappearance of strong magnetism in the high- 
pressure phase is a natural part of the story. Theo- 
retically, superconductivity mediated by residual spin 
fluctuations emerges at high pressure and low temper- 
atures. The conflicting experimental reports 14,15,17 on 
superconductivity in CaFe2As2 may be consistent with 
the strong sensitivity of T c on microscopic parameters 
expected from theory and with sample inhomogencitics 
as proposed in Ref. 17. 

The strong first-order volume collapse has only been 
observed in CaFe2As2 - so what is special about this 
compound? While we cannot give a definite answer at 
this point, the current status of both theory and exper- 
iment suggests that CaFe 2 As 2 (i) is a particularly soft 
material with a small c-axis lattice constant, (ii) displays 
a large electron-phonon coupling, and (iii) is located in 
close vicinity to a magnetic-non-magnetic transition. In- 
deed, first-principles calculations 59 have reported a gi- 
ant magneto-elastic coupling. In particular, the mag- 



nitude of the Fe moment has been found to be cou- 
pled to the c-axis lattice constant, with this effect being 
strongest in CaFe2As2 as compared to other 122 and also 
to 1111 compounds. Moreover, a soft lattice is known to 
be a crucial ingredient to a first-order volume-collapse 
scenario. 38,41-44 Note that, in our calculations, a small 
volume jump remains even for hard lattices due to the 
first-order nature of the purely electronic transition. 

Our conceptual ideas are not in conflict with itinerant 
spin-density-wave descriptions of the FeAs magnetism, 
but constitute a more strong-coupling-inspired view on 
the same physics. We speculate that an orbital-selective 
Mott scenario, likely with a continuous instead of a first- 
order transition, could apply to doping-driven transitions 
in iron arsenides as well, as a specific filling of the con- 
duction band is not required. 
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